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ABSTRACT 

Since the detection of very high energy (VHE) 7 -rays from Mrk 501, its broad 
band emission of radiation was mostly and quite effectively modeled using one 
zone emission scenario. However, broadband spectral and flux variability studies 
enabled by the multiwavelength campaigns carried out during the recent years 
have revealed rather complex behavior of Mrk 501. The observed emission from 
Mrk 501 could be due to a complex superposition of multiple emission zones. 
Moreover new evidences of detection of very hard intrinsic 7 -ray spectra obtained 
from Ferrm-LAT observations have challenged the theories about origin of VHE 
7 -rays. Our studies based on Ferrm-LAT data indicate the existence of two 
separate components in the spectrum, one for low energy 7 -rays and the other 
for high energy 7 -rays. Using multiwaveband data from several ground and space 
based instruments, in addition to HAGAR data, the spectral energy distribution 
of Mrk 501 is obtained for various flux states observed during 2011. In the present 
work, this observed broadband spectral energy distribution is reproduced with a 
leptonic, multi-zone Synchrotron Self-Compton model. 

Subject headings: BL Lacertae objects: individual (Mrk 501), Jets 
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Introduction 


The BL Lac source Mrk 501 (z=0.034) belongs to a sub-class of active galactic nuclei 
(AGN) that are known as high-energy peaked blazars (HBL). The broadband emission 
(radio to 7 -rays) of these objects is dominated by non-thermal radiation which is produced 
in the innermost part of the jets, that are oriented very close to our line of sight. This 
broadband emission is strongly Doppler-boosted. Like other TeV blazars, the spectral 
energy distribution (SED) of Mrk 501 characteristically shows a double-peaked profile. 
These peaks occur at keV and GeV/TeV energies when the SED is plotted in the vF„ 
versus v representation. The general understanding is that the first hump of the SED is 
caused by synchrotron radiation from the electron population gyrating in magnetic fields 
of the jet, but the origin of the GeV/TeV hump is unclear. The composition of these jets 
is also not known; it is not clear whether they are made of electron-positron plasma or 
electron-proton plasma. Even though Mrk 501 has been observed over last two decades in 
the entire electromagnetic spectrum, the existing multifrequency data could not provide 
explicit answer for the physical mechanisms that are responsible for the production of the 


GeV/TeV hump. T 


in lcptonic models (IKrawczvnski 


lis hump may be 


produced by interaction of e 


2004; 


Dermer fc Schlickeiser 


1993 


ectrons wit h ph otons 


Ghise 


1996 

) or protons with photon 

ields 

Miic 

te et al. 

2003; 

Mannheim 

1998) 


2012 1 ). 


lini fc Madau 


19981) or by mixed lepto-hadronic scenario O Cerrut i et al. 


,2000 


studied extensivelv in the past 

1999 

Krawczvnski et al. 

2000; 


by ( 

Sambruna e 

al. 

2000 ; 

Pian et al. 

1998; 

Villata & Raiteri 

Tavecchio et al. 

2001 

; Ghiscllini et al. 

2002 

), but the nature 


of this object is still far from being understood. The main reasons for this are the moderate 
sensitivities of the 7 -ray instruments, and the lack of simultaneous multifrequency data 
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during long periods. 


Since th e dete ction of this source at above 500 GeV by the Whipple observatory 


(i Quinn et al. 


199611 . it has shown several high states over the entire electromag netic 


spectrum and also few orphan TeV flares flAbdo et al. 


2011 


Neronov et al. 


20_12|). Mrk 501 


is also k nown for its m ajor , 


y-rays flCatanese et al. 


1997 


l ong timescale and short timescale 


Pian et al. 


1998 


Xue fe Cui 


2005 


lares in X-r ays an d VHE 


Albert et al. 


20070 . One of 


its historical o utbursts was observed in 1997 when the flux at energies above 1 TeV reached 


upto 10 Crab (lAharonian et al. 


1999b 


a|) • Foll owing this outburst, th e average flux of VHE 


y-rays dropped to 0.3 Crab during 1998-1999 (lA har oni an et al . 


200 l|). 


Rapid, and intra-night variab ility 


electromagnetic spectrum ( Gupta et al 


l as be en displ ayed 


2008 


by Mr k 501 over t 


Albert et al. 


2007! : 


Gupta et al 


re entire 


2012 ). 


Fast variability over time scales of minutes has been detected 


making the study of this object very interesting ( Albert et al 


during TeV orphan flares 


20071 ). There have been 


several mechanisms proposed for producing th e observed v ariability in the jet emission, 


rang ing from plasma me e 


Crusius-Waetze 


Begelman et al. 


ranis ms flKrishan fe Wiital 119941 b beamed radiation (e.g. 


V Lesch (119981) ). coherent instability in a compact emission region (e .g. 


(120081) ). misaligned minijets inside the main jet (e.g . 


(120 101)), jet deceleration (iGeorganopoulos & Kazanas 


2003 


Levi rson 


anisotropic electron beam directed alon g the iet (iGhisellini et al. 


blob inside the i et (blob-in-jet model: (iKatarzvriski et al 


Giannios et al. 


20071), wiggles in an 


20091) . relativistic plasma 


such as firehose (ISubramanian et al. 


2001 ])) and plasma instability 


20121 ) caused by anisotropic electron beam. 


Some fundamental questions regarding this source such as the content of its jet, 
location and mechanism of 7 -ray emission and the origin of observed variability are still 
not answered unambiguously. In an attempt to improve our understanding of the source, 
we present in this work a detailed study of multiwaveband data taken during 2011 January 
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to 2012 March using ground and space based instruments. We have fitted multiwaveband 
SED with multi-zone SSC model and discussed constraints on the physical parameters of 
Mrk 501. 


In this paper, we study the multiwavelength-multi-epoch behavior of Mrk 501 during 
the year 2011. Five differ ent fl ux states at different epochs along with a quiescent state SED 


observed by 


Abdoetah 


(20111 ) are modeled with a two zone SSC scenario and compared. 


Multi wavelength observations and analysis is presented in §2, and results are discussed 
in §3. A description of our new two zone model and modeling of six different flux states 
are described in §4. Finally we discuss implications of our two zone model on the blazar 
parameters considering Mrk 501 as an example in §5 . 


2. Multiwavelength observations and analysis 

VHE y-rays observations were made using High Altitude GAmma Ray (HAGAR) 
telescope array, at Hanle, India. In addition, archival data from Large Area Telescope 
(LAT) onboard Fermi , Proportional Counters Array (PCA) and All Sky Monitor (ASM) 
onboard RXTE, X-ray Telescope (XRT), Ultra-Violet/Optical Telescope (UVOT), and 
Burst Alert Telescope (BAT) onboard Swift , SPOL and Owens Valley Radio Observatory 


Tabic 1: HAGAR observations of Mrk 501 in 2010 and 2011 


Epoch 

Total duration 

Excess number of 

Mean y-ray rate 

Significance 


(min) (Pairs) 

ON source events 

(/min) 

a 

2010 March 22-2010 May 20 

400.2 (10) 

1577.0 ± 511.9 

3.9 ± 1.3 

3.1 

2011 March 31-2011 April 10 

279.5 (7) 

989.2 ± 399.6 

3.5 ± 1.4 

2.5 

2011 April 28-2011 May 10 

348.6 (9) 

2308.1 ± 454.9 

6.6 ± 1.3 

5.1 

2011 May 26-2011 June 03 

212.6 (6) 

967.1 ± 369.2 

4.6 ± 1.7 

2.6 
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(OVRO) were analyzed to obtain light curves and energy spectra. 


2.1. Optical and radio data 


The optical and radio data made available from Fermi multiwavelength support 
program are used. The optical observations were made by SPOL team using the SPOL 


CCD Imaging/Spectropolarimeter at Steward Observatory (jSmith et al . 20091). The optical 


V-band photometric and polarimetric fluxes are made publicly available by them on 
websiteqj and this data is used to obtain light curves and SEDs. 

The 15 GHz radio observations were made by using a 40 meter single-dish telescope 
at Owens Valley Radio Observatory (OVRO). The radio fluxes are also made publicly 
available by OVRO collaboration on websitedj and they are used t o obta in light cu rve at 


radio wavelengths. Details of the analysis method are described in (iRichards et al. 


coni)). 


In addition to this, we have also plotted radio data from Very Long Baseline Array (VLBA) 
at frequenci es 5 G Hz, 43 GHz and Submillimeter Array (SMA) at 230 GHz, obtained from 


Abdo et al. 


(20111) on each SED for the reference, as we do not have any radio observations 


of the core during 2011-2012. 


2.2. RXTE and Swift 


The PCA ( Bradt et al. 


19931 ) onboard RXTE is an array of five identical xenon-hlled 


proportional counter units (PCUs). The PCUs cover energy range from 2-60 keV with 
a total collecting area of 6500 cm 2 . We have analyzed standard 2 PCA data that have a 


Tittp://jarnes.as.arizona.edu/psmith/Fermi/ 
2 http: / / www. astro .caltech.edu / ovroblazars / 
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time resolution of 16 seconds with energy information in 128 channels. Data analysis was 
performed using HEASOFT (version 6.10). Data from PCA were analyzed to obtain the 
X-ray energy spectrum and light curve. For each of the observations data were filtered 
using the standard procedure provided in the RXTE Cook Book. The background models 
were generated with the tool “pcabackest”, based on RXTE GOF calibration hies for a 
’faint’ source (less than 40 ct/sec/PCU). 


The XRT onboard Swii 

"t uses a grazing incidence Wolter I telescope to focus X-rays 

onto a CCD ( 

Burrows et al. 

2005 

). The instrument has an effective area of 110 cm 2 , 23.6 


arcmin FOV, 15 arcsec resolution (half-power diameter), and covers an energy range of 
0.2-10 keV. The windowed timing (WT) mode data were used to obtain the spectrum 
from Swift-X RT. Source photons were extracted using a box region with the length of 40 


pixels and width about 20 pixels. Events with grades 0-2 were selected. The spectral data 
were rebinned by GRPPHA 3.0.0 with minimum 20 photons per bin. Standard auxiliary 
response hies and response matrices were used. 


A combined spectral ht was obtained for PCA and XRT data by normalizing the PCA 
spectrum with XRT spectrum. The PCA and XRT spectra in the energy range of 0.3-30 
keV were htted by using XSPEC with a cutoff powerlaw with linc-of-sight absorption. The 
li ne-of-sight absorptio n was hxed to a neutral hydrogen column density of 1.56 xlO 20 cm” 2 


flKalberla et al. 


2005!) 


Swift-XRT light curves are obtained from Fermi multiwavelength support program 


websiteqj and used in this study. 

The “Dwell” data from RATE-ASM were obtained from the ASM websitd 3 4 


and were 


3 http://www.swift.psu.edu/monitoring/ 

4 http://xte. mit.edu/ 













analyzed with the method discussed in IChi t nis et al . (120091) . A daily average flux between 
15-50 keV from Swift-BAT was obt ained from BAT websitjf], a detailed analysis procedure 


can be found in 


Krinnn ct al. 


(120131 


The Swift -UVOT (IRoming et al. 


20051) data were used to obtain fluxes in UVW1, 


UVM2 and UVW2 filters for different epochs. The snapshots of every individual observation 
were integrated with uvotimsum task and then analyzed with the uvotsource task. A source 
region of 10" radius was selected around the source, while the background was extracted 


from a circular region of V which is centered in a source-free region. The flux ob t ainec 


corrected for Galactic extinction of E( B _ V )=0.02 mag as given by 
each spectral band. 


Schlegel et al. 


(119981) in 


was 


2.3. Fermi- LAT 


The Fermz-LAT is a pair production telescope (lAtwood et al. 


2009 ) onboard the Fermi 


spacecraft. LAT covers the energy range from 20 MeV-300 GeV with a field of view > 2.5 
sr. The Fermi -LAT y-ray data of Mrk 501 over the period of 450 days (MJD: 55560- 56010) 
were obtained from websitqj. Data above 200 MeV were analyzed using the standard 
analysis procedure (ScienceTools-v9r31pl) provided by the Ferrm-LAT collaboration. 


A circular region of 10° radius “region of interest (ROI)” was chosen aroimd Mrk 501 
for event reconstruction from the so-called “diffuse” event class data which has the 
maximum probability of being the source photons. Events having a zenith angle <100° are 
only retained to avoid the background from Earth’s albedo. The spectral analysis of the 
resulting data set was carried out by including galactic diffuse emission component model 


5 http: //heasarc. nasa.gov/docs/swift / results/transients/ 
6 http: / / fermi.gsfc.nasa.gov/ 
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(gal_2yearp7v6_v0.fits) and an isotropic background component model (iso_p7v6source.txt) 
with post-launch instrumental r esponse functio n P7SOURCE_V6, using unbinned maximum 


likelihood analysis ((Cash 


1979 


Mattox et ah 


19961 ). A power law was used to model 


the source energy spectrum above 200 MeV, with integral flux and photon index as 
free parameters. The flux, and spectrum were determined by using unbinned GTLIKE 
algorithm. 


2.4. HAGAR 


The High Altitude GAmma Ray (HAGAR) telescope array is an Atmospheric 
Cherenkov Telescope array used for detecting VHE y-rays from the celestial sources. This 
array uses the wavefront sampling technique, and is located at the Indian Astronomical 
Observatory (IAO), Hanle (32° 46' 46" N, 78° 58' 35" E), in the Ladakh region of India, 
at an altitude of 4270 m. HAGAR consists of an array of seven telescopes arranged in 
the form of a hexagon, with one telescope at the center. Each telescope is separated by 
50 m di stance from its neighboring telescope. These telescopes use alt-azimuth mounting 


system (jGothe et al. 


20131 ). All seven telescopes have seven para-axially mounted front 


coated parabolic mirrors of diameter 0.9 m in each, with a UV-sensitive photo-tube at 
the focus of individual mirrors. These parabolic mirrors have f/D ratio of 1 and they 
were fabricated by using 10 mm thick float glass sheets. FOV of HAGAR telescope is 3° 
FWHM. The photomultiplier tubes (PMT) which are mounted at the focus of these mirrors 
are manufactured by Photonis (XP2268B) and have a peak quantum efficiency of 24% at 
400 nm. The high voltages fed to these PMT are monitored and controlled by C.A.E.N 
controller module (SY1527). In addition to recording signals from individual PMTs, the 
signals from the 7 PMTs of a telescope are linearly added to form a telescope pulse which 
is also recorded. The presence of any 4 telescope pulses above a preset threshold value 
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and within a window of about 150 ns forms the trigger for initiating data acquisition. The 
typical trigger rate was about 12 Hz. 


The CAMAC based Data Acquisition (DAQ) system is used in HAGAR. Relative 
arrival time of the Cherenkov shower front at each mirror is recorded for each event, as 
measured by TDCs with a resolution of 250 ps. The Cherenkov photon density at each 
telescope is measured by the total charge present in PMT pulses using 12 bit QDCs and 
a Real Time Clock (RTC) module synchronized with GPS is used to record the absolute 
arrival time of these events accurate upto /zs. In addition to this, a parallel DAQ using 
commercial waveform digitizers with a sampling rate of 1 GS/s (ACQIRIS make model 
DC271A) is also used to record telescope pulses. 


The energy threshold of the HAGAR telescope array is estimated to be 208 GeV 
for vertically incident q-ray showers for a > four-fold trigger condition, for which the 
corresponding collection area is 3.44xlO 8 cm 2 . The corresponding sensitivity is such that 


HAGAR wil 
observation 


detect a C rab-nebula-like source at a significance level of 5cr in 17 hours of 


Saha et al 


20131) . with no additional criteria for the rejection of background 


cosmic ray events. 


Observations of Mrk 501 were made during March-May 2010 and March-June 2011 on 
moon-less, clear nights using HAGAR telescope. Observation details are provided in Table 
H] The observations were carried out by tracking the source or background region with 
all seven telescopes. Each source (ON run) was followed (or preceded) by a background 
(OFF run) with the same exposure time (typically 40 minutes) covering the same zenith 
angle range as that of the source to ensure that observations were carried out at almost 


the same energy threshold. Selection criteria were applied to i dentify goo d qua 


Data were analyzed according to the procedure discussed in 


Shukla, et al. 


ity data. 


( 20121 ). In this 


procedure , Cherenkov shower front is approximated with plane front and space angle i.e. 
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angle between shower axis and pointing direction of telescope is estimated. 7 -ray signal is 
estimated comparing space angle distribution from ON-OFF pair. Only events with signals 
in at least five telescopes (> five-fold) were analyzed to reduce systematic errors, which 
corresponds to an energy threshold of 250 GeV for 7 -rays. 


2.5. ARGO-YBJ 


The ARGO-YBJ experiment situated at the Yangbajing Laboratory, Tibet at 4300 nr 
a.s.l. was designed to study cosmic 7 -radiation, at an energy threshold of ~ 100 GeV, by 
means of the detection of small size air showers. ARGO-YBJ consists of a single layer 
of Resistive Plate Counters (RPC) detectors covering an area of ~ 6700 m 2 to detect air 
showers. We have used published TeV y-r ay data o f M rk 501 from ARGO-YBJ, collected 


during 2011 October 17 to November 22, flBartoli et ah 


201211 


3. Results 

Mrk 501 was observed during year 2010 and 2011 when it was in a moderate state 
of activity using HAGAR telescope and VHE 7 -rays were detected from it. The source 

Table 2: Time periods used to obtain SEDs 


Epoch 

MJD 

Dates 

Days 

SI 

55651 - 55661 

2011 March 31 - 2011 April 10 

10 

S2 

55679 - 55692 

2011 April 28 - 2011 May 11 

13 

S3 

55707 - 55716 

2011 May 26 - 2011 June 03 

9 

S4 

55860 - 55890 

2011 October 26 - 2011 November 25 

30 

S5 

55919 - 55934 

2011 December 24 - 2012 January 08 

15 
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was detected with a 5 o significance during 2011 May (MJD: 55679-55692) when average 
flux reached a peak flux of ~1.5 Crab units (1 Crab unit—4.2 counts/min for at least five 
telescopes triggering). The average integral flux in this observation period above 250 GeV 
is found to be 4.04 xlO ~ 10 ph/cm 2 /sec. The light curve based on HAGAR observations 
during 2011 March-June is shown in Figure ED The TeV 7 -rays from Mrk 501 are detected 
at a total significance of 6.7cr over the observation period of two years. The source had 
brightened up moderately during 2011 over the entire electromagnetic spectrum. The 
multiwavelength light curves, from radio to 7 -rays are used to understand its flux levels and 
variability. The multiwaveband quasi-simultaneous light curve of Mrk 501 during 2011-2012 
is plotted in Figure [2j A few moderate and high states are identified during this period 
to study the spectral variation with the activity in X-ray and 7 -ray bands. SEDs were 
obtained for five such states by fitting multiwavelength data with SSC model. Details of 
these states are provided in Table [2] 


3.1. Flux and spectral variation during 2011-2012 

Multiwavelength flux and spectral variability of Mrk 501 as measured and reported 
by several ground and space based instruments during the 2011-2012 are presented in 
this section. The multiwaveband quasi-simultaneous light curve of Mrk 501 during 2010 
December 30 to 2012 March 24 is shown in Figure [2] The panels, in a descending order, 
correspond to data from OVRO (15GHz), SPOL (Optical V band), Swift-X RT (2-10 
keV), Swift -BAT (15-50 keV), Ferrm-LAT (0.2-2 GeV), and Fenro-LAT (2-300 GeV). 
The bottom panel corresponds to HAGAR data above 250 GeV. 

Mrk 501 was found to be variable in all the wavebands during the time span of 2011 
January to 2012 March with a few active states during this period. A clear variation of 
flux over a period of an year is observed in the radio, optical, X-rays, and 7 -rays. We 
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found that the source was brightest in X-rays and 7 -rays at the end of the year when it 
showed couple of flares in X-ray and 7 -rays. During this period of active states a few X-ray 
flares were observed by Swift -BAT (15-50 keV). Source showed flaring behaviour during 
2011 October 26-2011 November 25 and the peak flux in this duration was observed on 
MJD:55873, however peak flux observed by Swift-X RT was on MJD: 55931.09. The X-ray 
(Swift- BAT) light curve shows mild correlation with high energy 7 -ray light curve (2-300 
GeV). This correlation appe ars to be stro nger in high state of the source. A similar trend 


was reported earlier also by 


Gliozzi et ah 


(120061 ). Fluxes in different wavebands during 


HAGAR observation periods are also plotted, see Figure El Source was bright during month 
of 2011 May in entire electromagnetic spectrum. 


The Fermi -LAT light curves at different energies (0.2-300 GeV, 0.2-2 GeV, and 2-300 
GeV) during the period from 2010 December 30 to 2012 March 24 are plotted with a 
bin size of 15 days in Figure |4j The top two panels of this figure correspond to flux and 
photon index of low 7 -ray energy band ( 0 . 2-2 GeV), next two panels correspond to flux 
and photon index of high 7 -ray energy band (2-300 GeV). The bottom two panels of this 
figure correspond to flux and photon index of full 7 -ray energy band (0.2-300 GeV). Photon 
indices are only plotted for those bins that have more than 5<r detection. The LAT observed 
flux and spectral variation in low energy 7 -ray band ( 0 . 2-2 GeV) is found to be different 
compared to high energy band (2-300 GeV). 


Studies of spectral properties of Mrk 501 in X-ray and 7 -ray bands show significant 
spectral variability during the year 2011 in both the bands. Comparison of spectral variation 
with the activity of the source was carried out for the period (MJD: 55651-55934) using 
RXTE-FCA, Swift-X RT and Fermi- LAT instruments. X-ray and 7 -ray spectral indices 
were obtained for several flux states and given in Table [3l This table contains time interval, 
7 -ray flux above 0.2 GeV, 7 -ray photon index, and values of likelihood test statistics (TS) 
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as obtained from Fermi-LKY analysis for which spectra are made and last column contains 
X-ray indices as obtained from combined fit of Svnff-X RT and RXTE-PCA data. 


The Fermi -LAT collaboration had reported the value of photon index to be 1.78 based 
on the first 480 days of their observations of Mrk 501. This refers to an average spectrum 
mostly during the quiescent state. They detected remarkable spectral variability where the 
observed spe ctral index ranges from hardest value of 1.52±0.14 to the softest 2.51±0.20 


( Abdo et ah 


2 0111) . This change in the spectrum was not found to be correlated with the 


measured flux variations above 0.3 GeV. We also have detected large 


spectr a l var iabili 


during year 2011 in our study. We found much harder spectra than flAbdo et ah 


20111 ) 


y 


during our study, hardest spectrum was detected during S3 with the value of 1.27±0.21 in 
energy band of 0.2-300 GeV, (for details see Table [3]). Using 15 days bins, we have also 
detected very hard spectra of indcies ~ 1.2 during our 2011-2012 analysis, at higher energy 
y-rays band, (see fourth panel Figure |4l). The spectral variability is also seen in X-ray band 
during 2011. 


The cross plot between y-ray flux and photon index, shown in Figure 0 clearly shows 
two populations, one for low energy bins of 0.2-2 GeV and second for 2-300 GeV. The 
photon index increases with the increase in flux for lower energies (plotted in green square), 
but the cross plot shows scatter in case of high energy y-rays and no significant trend is 
visible (plotted in red downward triangles). This property indicates that low energy y-rays 
may be produced in a different emission zone having slightly different electron energy 
distribution and magnetic field than higher energy VHE y-rays in Mrk 501. 


In addition to the results discussed in this paper, the recent Ferm i -LAT observa tions of 

). New 


Mrk 501 challenge our present understanding about this source flNeronov et al. 


2012 |) 


evidence for the presence of very hard intrinsic y-ray spectra obtained from Fermz-LAT 
observations has challenged the theories of origin of VHE y-rays. Several very interesting 
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and viable explanations for the observed hard spectra have been proposed in the recent 


years. A hard 7 -ray spectrum could be obtained muc 


as discussed in proton synchrotron model (jAharonian 


1 more easily in hadronic scenario 


200Q|). Whereas, achieving hard 


spectrum from leptonic models is more demanding. Some viable scenarios, such as, 


relativistic Maxwellian-type electron energy dis 


acceleration process as cause of hard spectra flLefa et ah 


ributions which a re formed by a stochastic 


201 lb l) or the hard spectra 


produced by an electro magneti c cascade initiated by very-high-energy 7 -rays in the 


intergalactic medium flNeronov et ah 


2 0121 ) are available in the literature. 


4. Spectral Energy Distribution 


Since its detection at VHE 7 -rays, the broad-band emission from Mrk 501 has been 
explained mostly using one zone models in liter ature so far, which was quite effective to 


expla in the observed data from radio to 7 -rays (ISambruna et al 


2000 


Kra wcz vnski et al 


20Q0I). The main reason of the apparent success of one zone model was the lack of 


the simultaneous multiwaveband data and absence of data at low 7 -ray energies. The 
broad-band spectral and flux variability studies enabled by the multiwavelength campaigns 
carried out during the recent years revealed rather complex behavior of Mrk 501. The 
modeling of Mrk 501 SED assuming a homogeneous single emission zone is a very simplified 
situation. The broad band emission from blazar may be produced in an inhomogeneous 
region. The observed emission from Mrk 501 could be d ue to a co mplex superposition 


of multiple emission zones. The study carried out by 


Neronov et al. 


( 20 T 2 ) and the one 


presented here indicates the existence of two separate components in the spectrum, one 
for low energy 7 -rays and one for high energy 7 -rays. For example, emission from single 
zone is not enough to explain broad band emission of Mrk 501 as shown in Figure [ 6 ] where 
Mrk 501 SED is modeled using single zone. 
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We have considered two zone scenario to explain the broad band SED of Mrk 501. We 
assume that the observed broad-band SED is sum of two components (two zones) which are 
radiating simultaneously and boosted with almost same Doppler factor. These two zones 
are having comoving radii R* n and R oli f, and travel with a bulk Lorentz factor T towards 
the observer. The emission zones are filled with randomly oriented uniform magnetic fields 
B in and B out and isotropic population of non-thermal electrons. The energy spectra of the 
injected electrons in the jet frame are described by broken power laws with low-energy 
(E min to Efo) and high-energy (E& to B max ) components with indices of pi and p2. The 
outer zone is responsible for quiescent state flux and other compact inner zone which is 
close to the black hole is responsible for flaring activity in the jet. The broad band emission 
from blazar zone is sum of flux of quiescent component and active component. 

The radius of the emission zone is constrained by the variability time scales. The 
comoving radius of the emission zone is defined as 


R ~ cSt var / (1 + z). 


( 1 ) 


The values of t var for both zones are provided in Table [4l These values are consistent with 
th e flux variab ility observed in Mrk 501 during our study period and available in literature 


f flAcciari et al.l 201111 and references there in). 


Five different flux states have been identified during 2011, when the source was in 
moderate and high state of activity (see Figure [2] and Table [2]). An average TeV spectrum as 
observed by ARGO-YBJ during October 17 to November 22 is used to model S4 state. The 
SED obtained as a result of the 4.5-month-long multifrequency campaign (2009 March 15- 
2009 August 01 j o rganized by Fermi and MAGIC collaboration (here after MAGIC SED) 


flAbdo et all 


2 0111 ) is also modeled with two zone model for comparison. Multiwaveband 


data for all six SEDs (SI to S5 and MAGIC SED) are fitted with two zone SSC model 
and satisfactory fits are obtained, (see Figures [T1-[T2]). One more point to be noted is that 
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BAT flux covering the energy range of 15-50 keV during SI state is somewhat higher than 
the fitted model, whereas fitted model agrees well with the data from Swift-XRT and 
RXTE-PCA together covering the energy range of 0.3-30 keV. For other states, BAT flux 
seems to roughly agree with the fitted model and measurements from lower energy X-rays’. 

No change is observed in the jet flow (Doppler factor) during 2011 in fitted SED. 
Magnetic field strength in outer zone shows no variation among these states but, SED 
modeling of inner zone indicates that magnetic field in this zone varies with state. We found 
comparatively lower magnetic held strength during active states S2 and S4 in inner zone, 
but found higher magnetic held in high state S5 by SED modeling. We have also not seen 
any significant change in electron energy distribution of outer zone but change is indicated 
in inner zone. Higher electron energy density is found at the time of activity in inner zone. 
Also, we noticed significant change in electron energy spectral index after the break (p2) in 
this zone. 

Table 3: Fermi -LAT spectrum 


Epoch 

Flux (0.2 - 300 GeV) 

x 10 -8 ph/cm 2 /sec 

Photon Index 

0.2-300 Ge V 

TS 

XRT-PCA 

Index 

SI 

3.8T0.3 

1.59 ± 0.14 

102 

1.95 

S2 

5.12T0.88 

1.66 ± 0.13 

126 

1.87 

S3 

1.784=0.34 

1.27 ± 0.21 

67 

2.11 

S4 

4.5=1=0.23 

1.62 ± 0.09 

297 

1.75 

S5 

5.9=1=0.45 

1.46 ± 0.10 

276 

1.59 
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5. Discussion 


Mrk 501 is a core dominated radio source, with a one sided jet on pc scale, which 
extends till ~500 pc. This jet shows several sharp bends followe d by rap i d exp ansion and 


limb brightening structures on physical scale of ~lpc as seen by 


In recent years 


Giroletti et al. 


(120041 1 


s several attemp 

;s were made to mode 

the SED of Mrk 501 with 

(Ghisellini et al. 

2005; 

Graff et al. 

2008 

Giannios et al. 

2009j; 

Lefa et al. 


20111 ) 


2011 al) or adding extra break in the injected electron distribution ( Abdo et al. 
Alternatively, the SED could also originate from a two- component spine-sheath structure of 


the jet transverse to its direction 
radio jet morphology of Mrk 501 


Ghisellini et al 


Giroletti et al. 


20051) . as suggested by the complex VLBI 


20041 ). However, the energy transport in 


blazar jets generally occurs along the jet axis as witnessed by the motion of VLBI radio 
knots suggesting to logically connect the flux variability with the longitudinal evolution of 
components. 


In this section, we discuss the implications of the physical parameters of the source 
resulting from the SSC modeling of the SED. We also try to understand properties of 
the electron energy distribution emerging from SED modeling and constrain the physical 
processes responsible for the particle acceleration. We also examine the broadband 
variability of Mrk 501 in the framework of our two zone model. 


5.1. Variability 

An alternative way to constrain the physical parameters of the jet is to model its flux 
variability. Measurement of fast and rapid flux variability can shed light on motions of the 
bulk outflows of the plasma in innermost region of jets which is well beyond the current 
imaging capabilities of telescopes in any part of the electromagnetic spectrum. 
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In the SSC scenario the highest energy tail of the electron energy distribution (y>y br) 
is responsible for the production of the observed X-ray synchrotron continuum at > 0.5 
keV in HBLs, while the TeV 7 -rays might be produced through upscattering of synchrotron 
photons by the same population of electrons. The observed optical and X-ray variability 
during 2011-2012 may be explained by injection of fresh electrons in emission zones and 
cooling of the electrons due to SSC mechanism. Mrk 501 also shows energy dependent flux 
and spectral variability in 7 -rays. The source flux varies differently from lower energies 
(0.2-2 GeV) to higher energies (2-300 GeV). The observed 7 -ray variability is mainly 
divided into two bands, < 2 GeV and above 2 GeV. The 0.2-2 GeV 7 -rays observed 
by Fermi -LAT could be produced by low-energy electrons through IC scattering of UV 
synchrotron photons. The observed HE (>2 GeV) 7 -rays by Ferrm-LAT and VHE 7 -rays 
by HAGAR could be produced by IC scattering of the electrons having a Lorentz factor in 
the range ~ 10 4 — 10 5 . 

Mrk 501 was detected in moderate activity state during May 2011 by HAGAR as seen 
in last panels of Figures [2] and |3l A positive correlation between low energy X-rays and 
7 -rays is seen during this period, with the peak flux being observed in May 2011 in the 7 -ray 
and X-ray wavebands. Observed flux enhancement during the S2 period can be explained 
by the injection of fresh electrons in active zone (inner zone) of jet. These electrons may 
be accelerated to higher energies by shock in the jet. Two zone spectral modeling of this 
period suggest that the injected electron spectrum in inner zone has a power-law of index 
2, which could arise due to Fermi first order mechanism. Moreover, the Swift -BAT light 
curve is found to be anti-correlated with the other wavebands and we can not conclusively 
say why the flux in Swift -BAT band is found to be high during SI. 
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5.2. Application of two zone model in Mrk 501 multiwaveband data during 

2011 

The multifrequency data set provides an opportunity to obtain broad band SED of 
Mrk 501 in the quiescent and moderate states, and it also allows comparison between 


these flux states. In the work presented here we have compared five different 


of Mrk 501 during the year 2011, along with quiescent state S 


020111) . We have modified the single zone model developed by 


ID observed by 


lux states 


Abdoetah 


Krawczvnski et al 


(120041 ) to 


a multi zone model and used it to explain broad band emission from Mrk 501. Details of 


this one zone model can be found in 
by 


Shukla et al. 


Krawczvnsk i et al. 


(120041) . Same model was also used 


( 2012 1 to explain SED of Mrk 421. 


We have found that y-ray emitting zones are very close to the black hole around ~0.08 
pc which is consistent with other results presented in literature. The inferred magnetic 
field from the modeling of SED by our model is also in good agree ment with the m agnetic 


field claimed for the partially resolved radio core of Mrk 501 flGiroletti et al. 


200J). But 


m 


our model, source is not found in equipartition with the relativistic electrons along with 
magnetic field. On the other hand, source is consistent with being in equi part ition with the 


relativistic electrons in in spine-sheath model given by 


Ghisellini et al. 


(120051 ). From our 


study of multiwaveband SED modeling, we also infer that plasma in the jet is moving with 


Doppler factor of ~ 12 and it does not c hange w i 


consistent with other previous works (lAbdo et al. 


i the flux state and time. This value is 


20111 ) . Also we have not detected any 


change in the speed of the jet in y-ray emitting zone (<0.1 pc). The electrons in these blobs 
are accelerated through Fermi first order mechanism producing a power law distribution. 
We have found the inner blob has much narrower electron distribution than outer blob, 
with high minimum cutoff for y min . This blob is responsible for activity in the blazar zone 
and also the observed hard spectrum in the source. A narrow electron distribution can 
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produce very hare 


Katarzvriski et al 


spectr u m, similar s uggesti ons were also made by (ITavecchio et al.l (120091) ; 


(2006); 


Lefa et al 


(1201 lbl )). Electron population of outer zone is evolved 


and old, and this population has suffered radiative and adiabatic losses. We have observed 
via SED fitting that contribution of outer zone is not constant and it varies with flux 
state. Among the six flux states we discussed in this work, we found outer zone contributes 
significantly to y-ray hump in all the states except S5. This could be possible only if outer 
zone of S5 is not left with high energy electrons at the time of activity. We found that at 
the time of activity inner blob becomes dominant and it may produce hard spectrum at 
highest energies. The observed SED of Mrk 501 is sum of total radiation emitted by two 
zones and the shape of SED depends on the relative contributions from each zone. If both 
the zones contribute to the total observed flux then SED might be observed with a plateau 
at lower energy y-rays and a hard spectrum at highest energy. The difference between the 
spectral indices below and above the break energy Ap=p 2 -pi determined by SED modeling 
are close to 1 in the case of SI, S2, S3 and MAGIC SED for inner zone, see Table SI The 
value, Ap=l is expected as a result of classical synchrotron cooling break for a uniform 
emission region. 

The observed LAT spectrum during 2011 May when the source was in moderate bright 
state, shows plateau in the SED at lower energies ~(0.2-5 GeV), and a break in the slope 
at GeV (see Figure [9]). A hard spectrum with photon index of 1.5 is detected in 
higher energy (2-300 GeV) band during these observations. A similar behaviour was also 
rep orted d uring a fla re observed in the first 480 days of Ferrm-LAT operation (2008-2009) 


in (lAbdo et a l. 


2 0111) . Spectrum of this flare was found to be very hard/flat (~1.1) in the 


10-200 GeV range. The SED in the 0.3-200 GeV r ange during this flare also shows a break 


in the slope, around 10 GeV (INeronov et al. 


20121 ). Very hard spectra (<1.3) are detected 


in Ferrm-LAT data, few times during 2011-2012 in low energy (0.2-2 GeV) as well as in 
high energy band (2-300 GeV), (see Figure [4]). Origin of this very hard spectrum is still 






















23 


under debate. 

This work used results provided by the ASM /RXTE teams at MIT. This study also 
used Swift /BAT transient monitor results provided by the Swift /BAT team. This research 
has also made use of data obtained from the High Energy Astrophysics Science Archive 
Research Center (HEASARC), provided by NASAs Goddard Space Flight Center. Data 
from the Steward Observatory spectropolarimetric monitoring project used supported by 
Fermi Guest Investigator grants NNX08AW56G and NNX09AU10G. Radio data at 15 GHz 
is used from OVRO 40 M Telescope and this Fermi blazar monitoring program is supported 
by NASA under award NNX08AW31G, and by the NSF under award 0808050. This 
research has made use of the XRT Data Analysis Software (XRTDAS) developed under 
the responsibility of the ASI Science Data Center (ASDC), Italy. We are grateful to the 
engineering and technical staff of IIA and TIFR, who have taken part in the construction 
of the HAGAR telescopes and contributed to the setting-up of the front-end electronics 
and the data acquisition. We also thank Dr. David Paneque and Dr. Chen Songzhan for 
providing published data of Mrk 501 to use in this study. 


Tabic 4. SED parameters 


State 

T_var 

Magnetic 

Doppler 

1 l°g E min 

2 log E/max 

log E break 

Pi 

p2 

4 U e 

5 ^ 

- 

(hr) 

field 

factor 

[eV] 

[eV] 

[eV] 

- 

- 

[10- 1 2 3 4 5 ] 

[ u ’J u 'b\ 

- 

- 

(G) 

(-5) 

- 

- 

- 

- 

- 

(erg/cc) 

- 


MAGIC (Outer) 

48 

0.032 

12.07 

8.6 

11.6 

10.10 

2.4 

3.50 

1.8 

44.2 

SI (Outer) 

48 

0.028 

12.07 

8.9 

11.6 

10.10 

2.4 

3.95 

1.0 

32.0 

S2 (Outer) 

48 

0.028 

12.07 

8.9 

11.6 

10.10 

2.4 

3.95 

1.6 

51.3 

S3 (Outer) 

48 

0.028 

12.07 

8.9 

11.6 

10.10 

2.4 

3.95 

1.0 

32.0 

S4 (Outer) 

48 

0.028 

12.07 

8.9 

11.6 

10.10 

2.4 

3.95 

1.25 

40.1 

S5 (Outer) 

48 

0.028 

12.07 

8.9 

11.6 

10.10 

2.4 

3.95 

1.0 

32.0 

MAGIC (inner) 

6.9 

0.08 

12.07 

10.0 

11.95 

10.80 

2.0 

3.05 

11 

43.2 

SI (inner) 

6.9 

0.075 

12.07 

9.75 

11.85 

10.55 

2.0 

2.9 

18 

80.4 

S2 (inner) 

6.9 

0.056 

12.07 

10.2 

11.95 

10.45 

2.0 

2.85 

27 

216.0 

S3 (inner) 

6.9 

0.075 

12.07 

10.2 

12.00 

10.45 

2.0 

3.1 

20 

89.4 

S4 (inner) 

6.9 

0.056 

12.07 

9.7 

12.05 

10.90 

2.0 

2.5 

30 

240.0 

S5 (inner) 

6.9 

0.075 

12.07 

8.7 

11.90 

11.35 

2.0 

2.6 

32 

143.0 


1 E m i n : Minimum value of energy of the electrons present in the emission zone 

2 E max- Maximum value of energy of the electron present in the emission zone 

3 E break'- Break in the electron injection spectrum 

4 U e : Electron energy density 

5 77: Equipartition coefficient 
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Fig. 1.— HAGAR light curve of Mrk 501 during 2011 
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SI S2 S3 


S4 S5 



Fig. 2.— Multiwavelength light curves of Mrk 501 during 2011-2012. The bin size (in days) 
used for averaging the flux at different energy bands are mentioned at the right corner. The 
bin size for HAGAR observation are marked using x-error bars. The vertical dashed lines 
are representing the periods of moderate and high flux states for which SEDs are obtained, 
details of these states are provided in Table |2j 



































-27- 


1.05 
1.00 
13.86 
s=! 13.88 

— 13.90 
« 13.92 

— 13.94 
13.96 

6.5 

o 6.0 

CD 

(f) 

1 5.5 

g 

° 5.0 

4 '5 


0.02? f 
0.022 
0.020 
0.018 
0.016 
. 0.014 
0.012 
0.010 
0.003 


0.002 §■ 
0.001 \ 
0.000 §■ 


: 6 
; 4 
: 2 
b o 
, -2 
^ -4 
2.0 


® 1.5 r 


0.5 


S:8 

2.5 
2 . 01 - 

1.5 p- 
1.0 
0.5 p- 
0.01 

40 


OVRO (@ 15 GHz) 

i 




SPOL (Optical V Band* / Pol(%t 


t* ’ 


SWIFT-XRT (2 -10 keV) 


i 5 
=1 4 

M 

1 Q- 

==2 

_J 1 


RXTE-ASM (1.5 -10 keV) 



RXTE-PCA (2 -20 keV) 


= Swift-BAT (15-50 KeV) 


: Fermi-LAT (0.2-2 GeV) 


_ Fermi-LAT (2 - 300 GeV) 


HAGAR (> 250 GeV) 



80 100 
MJD (55600 +) 


Fig. 3.— Multiwavelength light curve of Mrk 501 during HAGAR observations 
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Fig. 4.— High energy 7 -ray light curve of Mrk 501 and photon index from Fermz-LAT during 


2011-2012 
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Fig. 5.— Cross plot: Flux vs Photon Index during 2011 with 15 days bin 
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Fig. 6.— One Zone SED of Mrk 501 during 2011 April-May [S2] 
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Fig. 7.— Two zone Spectral energy distribution of Mrk 501, as observed and presented in 
MAGIC SED. The solid black line show the best fits two zone model to the data, with the 
best-fit parameters listed in Table [4] 



Fig. 8.— Two zone SED of Mrk 501 during 2011 March-April [SI] 
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Fig. 9.— Two zone SED of Mrk 501 during 2011 April-May [S2] 



Fig. 10.— Two zone SED of Mrk 501 during 2011 June [S3] 




32 



Fig. 11.— Two zone SED of Mrk 501 during 2011 November [S4] 



Fig. 12.— Two zone SED of Mrk 501 during 2011 December [S5] 
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